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The effects of the existence of exotic nuclear shapes at the bottom of the neu- 
tron star inner crust - nuclear 'pasta' - on observational phenomena are estimated by 
comparing the limiting cases that those phases have a vanishing shear modulus and 
that they have the shear modulus of a crystalline solid . We estimate the effect on 
torsional crustal vibrations and on the maximum quadrupole ellipticity sustainable by 
the crust. The crust composition and transition densities are calculated consistently 
with the global properties, using a liquid drop model with a bulk nuclear equation of 
state (EoS) which allows a systematic variation of the nuclear symmetry energy. The 
symmetry energy J and its density dependence L at nuclear saturation density are 
the dominant nuclear inputs which determine the thickness of the crust, the range of 
densities at which pasta might appear, as well as global properties such as the radius 
and moment of inertia. We show the importance of calculating the global neutron star 
properties on the same footing as the crust EoS, and demonstrate that in the range of 
experimentally acceptable values of L, the pasta phase can alter the crust frequencies 
by up to a factor of three, exceeding the effects of superfluidity on the crust modes, 
and decrease the maximum quadrupole ellipticity sustainable by the crust by up to 
an order of magnitude. The signature of the pasta phases and the density dependence 
of the symmetry energy on the potential observables highlights the possibility of con- 
straining the EoS of dense, neutron-rich matter and the properties of the pasta phases 
using astrophysical observations. 
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oscillations - gamma rays: theory 



1 INTRODUCTION 

The neutron star crust is expected to be for the most 
part an elastic solid, comprising a Coulomb lattice of nu- 
clei, unbound electrons and, i n the inner crust, fluid neu- 
trons external to the nuclei (|Bavm et al.l Il971al lbl). How- 
ever, in the deep layers of the inner crust, near the crust- 
core boundary, the state of matter is more uncertain, and 
the possibility that the nuclei may undergo a series of 
shape transitions from spherical to cylindrical and slab, and 
then inverting to form the corres ponding bubble configu- 
rations, has been wid el y studied (iRavenhall et a l.l Il983bl: 



Hashimoto et aj 1984 : Ovamatsu et al. 1984: Ovamatsul 



1993*: 'Watanabe et al.1 l200ll : iMagierski fc HeenenI l2002f 
_Gogclcin ct al. 2 0071). More complex shapes can also exist 
l|Nakazato et al.ll2009l : iNewton fc Stonell2009l ') . These phases 
are collectively termed 'nuclear pasta'. Most studies have 
focussed on delineating the crustal region populated by the 
pasta phases. Important questions to address include: are 



there any practically observable neutron star phenomena 
that are sensitive to the properties of the pasta phase (as- 
suming they are different to those of a Coulomb crystal)? 
In what aspects of neutron star modeling do we need to se- 
riously pay attention to the properties of the pasta phases, 
and what modeling can we safely ignore the possibility of 
their existence? 

The sensitivity of the past a phases to neutrino i n- 
teractions have be en studied (|Horowitz et al.l l2004al lbl: 
,Gusakov et al.l |2004| ). and the possibility that the direct 
Urea cooling process is allowed in the bubble phases of the 
pasta introduces a potential observa tional probe should a 
young enough neutron star be found (|Gusakov et al.|[2004l ). 
The pasta phases are unlikely to be sol id; they are strongly 
analogous to terrestrial complex fluids (|Pethick fc PotekhinI 
1 19981 : i Watanabe fc Sonodall2005l ) which exhibit a wide range 
of responses to mechanical stimuli which are strongly de- 
pendent on the timescale over which the stimuli are applied. 
The modeling of mechanical properties such as the shear 
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modulus is a very complex task, and requires simulating the 
structure and mechanical properties of matter over a wide 
range of length scales and timescales (IPet hick fc PotekhinI 
Il998l ). Usually, in modeling phenomena associated with elas- 
tic properties of the neutron star crust, the shear modulus is 
taken to be that of a Coulomb crystal over the whole crust. 

In this paper we will estimate the maximal effect of the 
pasta phases on two elastic crustal processes that lead to 
potentially observable effects by comparing the case where 
the shear modulus in the pasta phases is taken to be that of 
a crystalline solid (which we refer to as 'solid' pasta), with 
the case where the shear modulus is set to zero (which we we 
refer to as 'liquid' pasta). The two neutron star properties 
we will examine are the maximal quadrupole deformation of 
the crust and the frequency of torsional oscillations of the 
crust. The former is important to determine the likelihood of 
observing gravitational wave e mission mechanism from ro- 
tating, deformed neutron stars (jAbbott et al.1120101 ') and the 
latter is a proposed mechanism for generating quasi-periodic 
oscillations (QPOs) in the tails of light cu rves of giant flares 
from soft gamma- ray repeaters (SGRs) ('Israel et al.' '2005); 
fWatts fc Strohniavcr 2006; Strohniayer_&JVatts 2005a, b). 

An important ingredient in the crust and core models is 
the symmetry energy of nuclear matter as a function of den- 
sity, i?sym(p), which can be defined as the energy difference 
per nucleon between pure neutron matter and symmetric nu- 
clear matter (equal numbers of protons and neutrons) at a 
given density. Its density dependence L at nuclear saturation 
density no ~ 0.16 fm~^, related to the pressure of pure neu- 
tron matter po through the relation L = po/3no, is strongly 
correlated with the crust-co re transition density and the 
thickness of the pasta layers (iHorowitz fc Piekarewicj l200ll : 
lOvamatsu fc lidal l2007l : | Xu et al.l |2008| '1. as well as global 
prope rties of the star such as the radius and moment of in- 
ertia l|Li fc Steineij [20061 : IWorlev et al.ll2"008l ). We will use a 
simple crustal model which is self-consistent with the core 
equation of state, and vary the value of the slope of the 
nuclear symmetry energy at saturation density over its ex- 
perimentally constrained range to take into account these 
correlations consistently. 

In section 2 we outline the approximate expressions we 
will use to estimate torsional frequencies and mountains, as 
well as give details about the EoS used and crustal model, 
including the shear modulus. In section 3 we give the results 
and in section 4 we discuss our conclusions. 



not confined to the nuclei. This is the zero-temperature ex- 
pression, a good approximation for neutron star tempera- 
tures below ~ 10*K. We assume that the shear modulus 
is isotropic. From this one can calculate the speed of shear 

waves in the crust «s = (/^/p) 2 , where p is the mass density 
corresponding to the baryon number density Wb. 

We will use the expression to calculate the following 
quantities: 

(i) The fundamental frequency and overtones of 
crustal torsional oscillations can be estimated as 
(|Samuelsson fc Andersson|[2007l ) 
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respectively, assuming an isotropic crust, n is the number 
of radial nodes the mode has, and I is the angular 'quan- 
tum' number. In order to estimate the superfluid effects of 
the dripped neutrons, one introduc es a multiplicative factor 
l|Andersson et al.ll2009l l. where 



(1-Xn 



[1 - X,l(m•/n^n)] 



(4) 



encodes the effect of entrainment of the free neutrons by 
the lattice, e^" = {l-2GM/c^R) and Rc is the radius of the 
crust-core boundary. This formula is derived from a plane 
wave analysis of the crustal shear perturb ation equations 
l|Pirol I2OO5I : ISamuelsson fc Andersso"nll200'it l. and has, in a 
related form and without superfluid effects, been used to 
study the i nfluence of the nuclear EoS on torsional mode 
frequencies (|Steiner fc Wattj|2009l ). The mesoscopic effec- 
tive mass m* encodes the 'drag' on the superfluid neutrons 
by the crystal lattice through which they fl ow; it has bee n 
calculated for a limited number of densities (|Chamelll2005l ). 
In this work we take the two extreme values m*/mn = 15 
and m^/m„ — 1, the latter corresponding to no entrainment 
effect. Note that the effective mass is the only quantity we do 
not calculate self-consistently with the crustal and core EoS 
used, something that should be addressed in future studies. 

(ii) An estimate of the maximum quadrupole ellipticity 
sustainable by the crust of a neutron star of mass M and 
radius R 



2 METHOD 

The shear modulus of a Coulomb lattice of positively 
charged nuclei in a uniform negatively charged back- 
ground in the neutron star crust at a baryon number 
density rtb was d etermined thro ugh Monte-Carlo simu- 
lation JOgata fc Ichim aru 199(|; IStrohmaver et aTl Il99ll : 
IChugunov fc Horowitz! [2Oi0i ) and can be written as 

^ ^ 0.1106 (^f^ A-*/^nt^^l - XO'/'(^e)^ (1) 

where the nuclei are characterized by the nucleon and 
proton number A, Z and X-a is the fraction of neutrons 



P {GM\ ^ _ 7 (o^^^\ 

where Vc is the crustal volume. This has the simple phys- 
ical interpretation of the ratio the stress energy to the 
gravitational energy of the star multiplied by the break- 
ing strain of the crust amax. This s imple estimate has been 
confirmed by rigoro us calculations (|Ushomirskv et al]|2000l : 
iHaskell et al.|[200^ ). The breaking strain has recently been 
estimated to be ~ 0.1 using molecular dynamics simulations 
(jHorowitz fc Kandaull2009h . Additionally, we calculate the 
equivalent gravitational wave strain amplitude from a ro- 
tating neutron star with a moment of inertia I, frequency u 
and distance to Earth d l| Abbott et al.ll2007l '): 
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From our crustal model, we need the crustal composi- 
tion and density at the crust-core and spherical-pasta tran- 
sitions. We us e the compressible liquid drop m odel (CLDM) 
developed in ()Newton. Gearheart fc LilbOllI '). Like the liq- 
uid drop and droplet models of terrestrial nuclei, the CLDM 
focusses on the average, macroscopic properties of nuclear 
clusters in the crust such as their mass and size, whilst ne- 
glecting quantum shell effects. The interior nuclear density 
is a free variable (i.e. the nuclear clusters are compressible) 
rather than specified a prion as the saturation density of 
symmetric nuclear matter. An external gas of dripped neu- 
trons, expected to coexist with the nuclear clusters in the 
inner crust, is calculated using the same bulk Hamiltonian 
as the matter inside the clusters. The interfacial energy be- 
tween the two bulk phases is obtained by specifying the 
magnitude and proton fraction dependenc e of the surface 
tension using the form of L attimer et al (iRavenhall et al.l 
Il983al : iLattimer et"al]ll985l ). The neutron star inner crust 
is assumed to have a periodic structure microscopically, and 
the Wigner-Seitz (W-S) approximation is employed in which 
the unit cell, which generally will have a cubic or more com- 
plicated structure, can be replaced by a cell of the same 
volume and possessing the same geometry as the nuclear 
cluster. We consider spherical, cylindrical and planar geome- 
tries for both normal nuclear clusters and bubbles. The equi- 
librium cell size, nuclear size, proton fraction and dripped 
neutron density are calculated by minimizing the total en- 
ergy of one cell at a given baryon number density. The en- 
ergetically preferred phase is determined by calculating the 
equilibrium energy density for each of the six nuclear geome- 
tries and for uniform matter, and selecting the lowest. By 
repeating the calculation for the density range of the inner 
crust, we obtain the composition and the transition densities 
(the transition to pasta is taken to be the transition between 
spherical-cylindrical geometries) . 

The bul k nuclear energy is given by a Skyrme-like inter- 
action MSL (|Chen et al.ll2009^ . The MSL interaction con- 
tains parameters directly related to the slope of the sym- 
metry energy, allowing L to be freely adjusted and giving 
us the ability to scan the appropriate parameter space of 
the nuclear interaction without changing the properties of 
the symmetric nuclear matter EoS. We take as a conser- 
vative raB£e_ofLftwnnucl«^ 25 < L < 115 



MeV (IChen et all 



20051; iLi fc Chen 2005; Shcttv et al.ll2007l: 



jgOOj ": 'Danielcwicz & Lcc 20'07'; iLi et al 
2009; Tsang et al. 2009; C hen et al 



Klimkiewicz et al 
2008; Cent cUes et al . 

2010l : IXu et al.ll2010l '). alt hough it is worth not i ng that the 
most recent constra ints l|Centelles et al. I l2009l : [Ch en et al.l 
I2OIOI : IXu et al.|[2O10l ) , coupled with infe rence from theoreti- 
cal calculations of pu re neutron matter (|Hebeler et al.|[2010l : 
iGandolfi et al.|[201ll ). place the value of L in the lower half 
of this range 25 - 75 MeV. We hold the absolute value of the 
nuclear symmetry energy at saturation density constant at 
J = Esyu^{no) = 32 MeV throughout. 

We use the same uniform nuclear matter EoS as input 



into the Tolman-Oppenheimer-Volkoff equations to calcu- 
late the global neutron star properties needed (mass M, ra- 
dius R, and crust thickness A = i? — 7?c, crust volume 14), 
and approximate the mo ment of inertia using the formula 
(jLattinier fc Schut jlioosi ) 
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As an upper estimate for the effect of nuclear pasta on 
the crustal shear phenomena, we will set the pasta shear 
modulus to zero. We will thus be exploring two extreme cases 
distinguished by the density nt at which we evaluate the 
shear modulus for use in Eqs. [2] and [5] 

(i) Pasta as an elastic solid; then the shear modulus is 
evaluated at the crust-core boundary nt = ncrust->corc 

(ii) Pasta as a liquid; then the shear modulus is evaluated 
at the boundary between the phase of spherical nuclei and 
the pasta phases nt = Wsph^pasta 

Taking pasta phases to be liquid, we are effectively mov- 
ing the boundary between the solid crust and the liquid 
core of the star down in density, so that the solid crust gets 
thinner. One would then expect the fundamental mode fre- 
quency to decrease, along with e. 



3 RESULTS 

Fig. 1 shows the density and pressure at the crust-core tran- 
sition and the transition from spherical nuclei to the pasta 
phases. The difference between the two transitions indicate 
the range of density and pressure inhabited by the pasta 
phases, versus L. The well-known inverse correlation be- 
tween the crust-core transition density and L is manifest; 
the corresponding inverse correlation between the crust-core 
transition pressure and L is also shown, although one should 
note th at this relation is dependent on the nuclear model 
used (iPucoin et al.|[201ll ') and is sensitive to other parame- 
ters of the nuclear EoS. The density range of pasta is largest 
at the smallest values of L, and disappears at the high end 
of the experimentally constrained range of L. 

Fig. 2 shows the crust composition parameters nec- 
essary for the determination of the shear modulus in a 
Coulomb crystal: the total nucleon number of the nuclear 
cluster A, the corresponding proton number Z and the frac- 
tion of dripped neutrons in the unit cell X-n. A and Z at 
the crust-core transition generally decrease as L increases, a 
trend that is largely due to the fact that the crust-core tran- 
sition density also decreases with L (and A and Z rise with 
density at the highest densities in the inner crust). At the 
spherical nuclei - pasta transition, A is reasonably uniform, 
while Z still decreases from 30 to about 10. It is impor- 
tant to note that in the semi-classical liquid drop model, 
shell effects are not included, which could significantly af- 
fect the trend. The dripped neutron fraction is larger 
at the spherical nuclei - pasta tran sition (as is the dr ipped 
neutron mesoscopic effective mass l|Carter et al.|[2005l )): Xn 
increases with L at both transition regions. 

Fig. 3 shows the shear modulus relative to the pres- 
sure, and the shear speed, at the two transition densities as 
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Figure 1. The transition density (top) and pressure (bottom) as 
a function of the slope of the nuclear symmetry energy L. Filled 
circles refer to the crust-core transition and empty circles show 
the transition between spherical nuclei and the pasta phases in 
the inner crust. 



Figure 3. The shear modulus compared with the pressure (top) 
and the shear speed (bottom) at the two different transition den- 
sities referred to in Fig. 1, as a function of the slope of the nuclear 
symmetry energy L. 
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Figure 2. Composition of the bottom of the inner crust: the 
total nucleon number (top), proton number (middle) and global 
dripped neutron fraction (bottom) at the two different transition 
densities referred to in Fig. 1, as a function of the slope of the 
nuclear symmetry energy L. 



a function of L, using the compositional and transition pa- 
rameters from figs. 1 and 2. It is immediately obvious that 
the influence of the pasta phases disappears at L > 80 MeV. 
The shear modulus and speed decrease as L increases; thus 
shear phenomena in the crust will be strongest for L at the 
lower end of the experimentally constrained range. 

Figs 4 and 5 show the relevant global neutron star and 
neutron star crust properties for a IAMq star. The radius 
R increases as L increases, as one would expect: since L at 
nuclear saturation density is a good fiducial indicator of the 
internal neutron star pressure, one would expect that the 
higher the pressure, the larger th e radius. Thi s is a relation 
that is well known (|Lattimer fc Pra kash 2 00ll : IlI fc Steineil 
l2006f ). The thickness and volume of the crust that contains 
spherical nuclei and the total crust with the pasta phases 
included both increase with L at L < 80 MeV and then 
decrease at higher L. This trend is the convolution of the 
decrease of the transition densities with L and the increase 
of the stellar radius with L. The moment of inertia follows 
a similar trend. The thickness of the pasta layer is greatest 
at intermediate values of L « 70 MeV, rising to about 10% 
of the crustal thickness and volume. It is important to note 
that, although the pasta layers occupy a small layer of the 
crust in terms of radius of volume, they occupy the highest 
density region of the crust, so the fraction of the mass of 
the crust occupied by the pasta layer s will be much higher, 
approaching 50% (|Lorenz et al.lll99^ 'l. 

Given the microscopic and global neutron star prop- 
erties, both consistently calculated with the same under- 
lying neutron star EoS, we can now estimate the maxi- 
mum quadrupole ellipticity sustainable by the crust and 
the fundamental mode frequency of torsional crust oscil- 
lations. Fig. 6 shows the maximum quadrupole ellipticity 
of the crust of a IAMq star, normalized to the canonical 
value of 10~^(a"brGak/10~^), as a function of L with the pasta 
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Figure 4. Global neutron star properties as a function of the 
slope of the nuclear symmetry energy L: for a IAMq neutron star, 
the radius (top), and moment of inertia (bottom) are plotted. 
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Figure 5. Global crust properties as a function of the slope of 
the nuclear symmetry energy L: for a IAMq neutron star, the 
crust thickness (top) and crust volume (bottom) are plotted. The 
latter two quantities are shown taking the bottom of the crust to 
be the true crust-core boundary (filled circles) and the spherical 
nuclei- pasta boundary (empty circles). 

phases taken to be a Coulomb solid and a liquid. The maxi- 
mum quadrupole ellipticity peaks at L « 50 MeV; the peak 
occurs because of the competing trends of increasing stel- 
lar radius and decreasing crust-core transition density as L 
increases, a feature which only emerges when the crust and 
core are treated consistently. For L < 70 MeV, and taking 
the pasta phases to be solid, the maximum ellipticity is close 
to the canonical value. However, if the pasta phases are as- 
sumed to be liquid, the maximum ellipticity is an order of 
magnitude lower ~ 10~*(a"brcak/10~^). At higher values of 
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Figure 6. Estimate of the maximum quadrupole ellipticity of 
a IAMq neutron star (top) and the corresponding gravitational 
wave strain (bottom) as a function of the slope of the nuclear 
symmetry energy L. The strain is calculated assuming a neutron 
star frequency of 300Hz and a distance of O.lkpc from Earth. The 
filled circles indicate the value taking the shear modulus in the 
pasta phases to be that of a solid Coulomb lattice; the empty 
circles indicate the value taking the shear modulus in the pasta 
phase to be zero (i.e. liquid pasta). The dashed line in the bottom 
plot indicates the se nsitivity of the most recent LIGO science run 
l|Abbott et al.ll2010l 'l. 

L, the pasta phases become insignificant, but the stiffening 
of the EoS and reduction of the crust-core transition density 
means that the crust thickness relative to the size of the star 
overall becomes smaller, and the maximum ellipticity sus- 
tainable falls by almost two orders of magnitude compared 
to the canonical value. The equivalent gravitational wave 
strain is plotted in the bottom half of Fig. 6 for a reference 
neutron star frequency of 300Hz and distance of O.lkpc. The 
sensitivity of the most recent LIGO science run to such grav- 
itational waves is indicated by the dashed line; this shows 
that, despite the reduced strain caused by a stiff EoS or the 
existence of liquid pasta phases, such GWs are still poten- 
tially observable with current or future detectors, offering 
the chance to distinguish between the various possibilities 
of crust properties. 

Fig. 7 shows the fundamental frequency I = 2,n — 
and first overtone / = 2, n = 1 of the crustal torsional os- 
cillation spectrum for a IAAIq star, where n is the num- 
ber of radial nodes and I the angular constant associated 
with the spherical harmonics Yi"^. Observed values of the fre- 
quency of QPOs at 16, 26, 28 and 30 Hz jls racl ct al. 200^; 
I Watts fc Strohmaveil [200^ : IStrohmaver fcWatts, ,2005a . bl). 
within the range of the calculated fundamental frequencies, 
are indicated by the dashed lines, as well as the observed 
84, 92, 150, 155 and 625 Hz frequencies within the range 
of the first overtone. The frequency generically decreases 
with increasing L as the crust becomes thinner relative to 
the radius of the star. Ignoring superfiuid effects and the 
effects of the pasta phases, the frequency matches observed 
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Figure 7. (Color online) The frequency of the fundamental tor- 
sional oscillation mode (top) and the first overtone (bottom) in 
the the crust for a 1.4Mq as a function of the slope of the nuclear 
symmetry energy L. The circles and triangles show the frequency 
assuming the shear modulus in the pasta phases to be that of a 
Coulomb lattice; the triangles take into account the entrainment 
of the supcrfiuid neutrons by the nuclear clusters. The squares 
and diamonds show the frequency assuming the shear modulus 
in the pasta phase to be zero (i.e. liquid pasta); the diamonds 
take into account the entrainment of the superfluid neutrons by 
the nuclear clusters. The dashed lines show possible candidate 
frequencies for the fundamental modes; 18, 26, 28, 30Hz in the 
fundamental frequency range and 84, 92, 150, 155Hz in the range 
of the first overtone. 



QPOs from SGRs only at the lowest values of L. If the pasta 
phases are liquid, the frequency falls by a factor of 3, making 
it difficult to match the 28Hz frequency observed, and be- 
ing consistent with the 18Hz observed frequency only at the 
lowest value of L. Compare to the effects of superfluidity, 
the effects of Itqmd pasta are about a factor of 4 larger. It 
should be noted that the observe d frequencies could also be 
explained by core Alfven waves (jSotani et ahllioosh . n = 
modes scale with [{I + 2)(Z — 1)]'^''^, so higher-? modes will 
show a similar effect of the pasta phases. The same trend is 
present for the first overtone. It is interesting to note that the 
observed 625Hz mode, often cited as the main candidate for 
the first overtone, is consistent only with low L < 50 MeV 
and solid pasta, whereas the lower observed frequencies 84- 
155Hz, which are often matched with n = Q, I > 2 modes, 
are consistent with a wider range of L, 60-110 MeV, and 
pasta properties. 



4 CONCLUSIONS 

We have demonstrated that the pasta phases can have a 
very significant effect on observable neutron star phenom- 
ena, reducing the frequencies of crustal torsional modes by 
up to a factor of 3 and the maximum quadrupole ellipticity 
sustainable by the crust by an order of magnitude. The ef- 



fect of the pasta phases is comparable with other transport 
properties such as entrainment of superfluid neutrons. In ad- 
dition, we have demonstrated that a consistent treatment of 
the crust composition and core equation of state is required 
in modeling these phenomena, and that when one uses such 
a consistent treatment the possibility of constraining the nu- 
clear symmetry energy at densities around saturation using 
astrophysical observations emerges. The models of torsional 
modes and mountain building in neutron star crusts have 
many other uncertainties. Having a diverse range of neu- 
tron star phenomena amenable to independent observation 
allows one to check the consistency of our models; study- 
ing the signature of nuclear matter properties in such phe- 
nomena opens up another set of constraints from nuclear 
experiment. 

Clearly, more work needs to be done investigating the 
likely structure of the pasta phases over a variety of length 
scales and estimating their actual shear modulus. This is 
likely to depend on the particular shapes present. One 
also cannot ignore temperature effects, especially for the 
long timescale proce ss of accreting material onto the crust 
(jChugunov fc Horowitz ,2010 ). The pasta phases are likely 
to be disordered on len gth scales of ~ 100 ce ll lengths due 
to thermal ffuctuations (jWatanabe et al.]|2000l ): however, we 
cannot rule out some ordering mechanism involving, for ex- 
ample, the magnetic fleld. One might envisage the cylinders 
or slabs aligning along a particular direction, giving the 
phase s liquid crystal-like properties (|Pethick fc PotekhinI 
Il998l ). The approximation of pasta being a liquid may be 
quite good for accretion processes which build mountains, 
as it is known that many complex fluids exhibit continuous 
flow in response to stress applied slowly over long timescales. 
A more microscopic treatment of crustal matter, including 
shell effects of nuclei and dripped neutrons, is also necessary 
for a more realistic determination of the transition densities 
and superffuid effects. We hope that demonstrating the fea- 
sibility of setting limits on those properties observationally 
provides additional motivation for these further studies. 

While finalizing this manuscript, we learned of a simi- 
lar study wh ich focusses on the torsional crust oscillations 
(|Sotanill201ll ). We briefly comment on how our study com- 
pares with the results obtained there. The main difference 
is that, whereas we use approximate methods for the de- 
scription of the torsional mo des and a co nsistent treat- 
ment of crust and core EoS, ISotanil (|201ll ) uses a more 
sophisticated calculation of the crustal frequencies, but a 
specific crust and core EoS which are not c onsist ent with 
each other. The core EoS used in ISotanil (|201ll ) is rela- 
tively stiff, giving ne utron star radii con sistent in our model 
with L > 100 MeV jSotani et al.ll2008l ). putting it outside 
of current best exp erimental and observational estimates 
(jSteiner et al.|[2O10l ). The crustal composition model is that 
of Douchin and Haensel (DH) IDouchin & Hacnsol 2001j), 
using the SLy4 EoS with L = 46 MeV. The transition den- 
sity to the pasta phases is taken as a free parameter rather 
than emerging consistently from the crustal model (indeed, 
the DH EoS predicts no pasta phases); the lower value taken, 
Tit ~ 0.007 fm~^, is well below that which is realistically ex- 
pected from liquid drop and more microscopic models. The 
range we find is 0.05 - 0.07 fm"^ = 7 x 10^^ - lO" g cm'^ 
The differences in our methods ma kes a d i rect c omparison 
of our results difficult. For example. ISotanil (|201ll ) obtains a 
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fundamental mode frequency of « 22 Hz for a IAMq neu- 
tron star when the shear modulus is taken to be that of a 
Coulomb crystal up to the crust-core boundary, compared 
with ~ 34 Hz in our model, taken at L = 46 MeV; however, 
the stiffness of the crust in the Sotani model would lead to 
a thinner crust and a decrease in the frequency compared 
with our results. It would be useful to make a comparison 
of the methods used using consistent nuclear physics. As a 
final note, we have estimated the comparative effect of the 
pasta phases and the effects of superfiu id neut r ons, f inding 
them to be comparable, as suggested in ISotanil (|201ll ). 
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